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ABSTRACT 

Broad absorption lines (BALs) in quasar spectra are prominent signatures of high- 
velocity outflows, which might be present in all quasars and could be a major contrib- 
(— i \ utor to feedback to galaxy evolution. Studying the variability in these BALs allows us 

O ■ to further our understanding of the structure, evolution, and basic physical properties 

of the outflows. This is the third paper in a series on a monitoring programme of 24 
luminous BAL quasars at redshifts 1.2 < z < 2.9. We focus here on the time-scales of 
variability in C iv A1549 BALs in our full multi-epoch sample, which covers time-scales 
from 0.02—8.7 yr in the quasar rest-frame. Our sample contains up to 13 epochs of 
data per quasar, with an average of 7 epochs per quasar. We find that both the inci- 
dence and the amplitude of variability are greater across longer time-scales. Part of our 
monitoring programme specifically targeted half of these BAL quasars at rest-frame 
time-scales ^2 months. This revealed variability down to the shortest time-scales we 
probe (8—10 days). Observed variations in only portions of BAL troughs or in lines 
that are optically thick suggest that at least some of these changes are caused by 
clouds (or some type of outflow substructures) moving across our lines of sight. In this 
crossing cloud scenario, the variability times constrain both the crossing speeds and 
the absorber locations. Specific results also depend on the emission and absorption 
geometries. We consider a range of geometries and use Keplerian rotational speeds 
to derive a general relationship between the variability times, crossing speeds, and 
outflow locations. Typical variability times of order ~1 year indicate crossing speeds 
of a few thousand km/s and radial distances near ~1 pc from the central black hole. 
However, the most rapid BAL changes occurring in 8-10 days require crossing speeds 
of 17 000 - 84 000 km s^ 1 and radial distances of only 0.001-0.02 pc. These speeds 
are similar to or greater than the observed radial outflow speeds, and the inferred 
locations are within the nominal radius of the broad emission line region. 

Key words: galaxies: active - quasars: general - quasars:absorption lines. 



1 INTRODUCTION 

High- velocity outflows are an integral part of the quasar sys- 
tem and likely contribute to feedback to the host galaxy. The 
material accreting onto the central supermassive black hole 
(SMBH) in quasar systems might release its angular momen- 
tum through these outflows. Moreover, these outflows might 
inject sufficient kinetic energy into the host galaxy to affect 
star formation, to aid in the 'unveiling' of dust-enshrouded 
quasars, and to aide in the distribution of metal-rich gas to 
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the in tergalactic medium (iDi Matteo et al1 [2005: Mol l et all 
l2007h . 

Many properties of these outflows, including their loca- 
tion and three-dimensional structure, are still poorly un- 
derstood. Current models of these outflows predict that 
they originate from the rotating accretion disc within the 
quasar system and are accelerated to hi gh speeds via radia- 



tive and /or magne t o-centrifugal force s dMurrav et all 



2005 



Proga et all l2000l : iProga fc Kallmanl 120041 : lEverett 
Progal 1200 7J). In order to test these models and estimate 
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the mass-loss rates and kinetic energy yields of these out- 
flows, improved observational constraints are necessary. Un- 
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derstanding these properties will help to illuminate the role 
of quasar outflows in feedback to the host galaxy. 

The most prominent signatures of these accretion disk 
outflows that are detected in quasar spectra are broad ab- 
sorption lines (BALs) . A quasar is classified as a BAL quasar 
if it has Civ Af549 absorption that reaches depths of least 
10 per cent below the continuum over >2000 km s _1 . 

This work is the third paper in a series on the vari- 
ability in these BALs. We introduce our ongoing mon i torin 
programme of 24 BAL quasars in ICapellupo et al.l (|201 
hereafter Paper 1). We begin with a sample o f quasars from 
the BAL variability study in IBarlowl l|l993h . and we have 
re-observed these quasars to obtain both a longer time base- 
line between observations and multiple epochs of data per 
object. In Paper 1, we directly compare Civ BAL variability 
properties between a "short-term" interval of 0.35—0.75 y 
and a "long-term" interval of 3.8—7.7 yr. In particular, we 
report that 65 per cent of quasars exhibited BAL variability 
in th e long-term, compar ed to 39 per cent in the short-term 
data. iGibson et al.l (|2008T ) detected variability in 92 per cent 
(12/13) of quasars on similar time-scales as our "long-term" 
data. This indicates that a majority of BAL quasars vary 
on multi-year time-scales. We also find a slight trend to- 
wards a larger change in strength in the long-term. These 
results are broadly consistent with the results of p r eviou s 
wor k on Civ BAL vari a bility , i.e. iLundgren et alj (I2007T ) 
and IGibson et all d200g|; l2010t). The second paper in our 
BAL variability series, ICapellupo et al l (|20ll hereafter Pa- 
per 2), includes all of the data obtained through 2009 March 
on our sample of 24 BAL quasars. We find that by adding 
more epochs of data, even at similar time-scales as the mea- 
surements in Paper 1, the fraction of quasars with variable 
BAL absorption increases. 

Paper 2 also discusses in detail the two prevailing sce- 
narios that are likely causing the observed BAL variability, 
i.e., a change in ionization in the outflowing gas and outflow 
clouds, or substructures, moving across our line of sight. 
Using the results of Paper 2, along with the results of pre- 
vious work, we can not rule out either scenario. Instead, 
Paper 2 concludes that the actual situation is likely a com- 
plex amalgam of both scenarios. Understanding the cause(s) 
of variability is necessary for using variability to constrain 
properties of the flows. In the moving-cloud scenario, the 
time-scales of variability can help constrain the sizes, speeds, 
and locations of clouds. For example, in Paper 1 we estimate 
that the distance of the outflowing gas is roughly within 6 pc 
of the central SMBH for those BALs that varied on time- 
scales of ~0.5 yr. Variability on shorter time-scales would 
place the outflow at even smaller distances, based on nom- 
inally shorte r crossing times for clouds moving across our 
line of sight (|Hamann et al.ll2008l ; Paper 1). In the changing 
ionization scenario, variability on shorter time-scales require 
high densities due to the shorter recombination times, which 
would also indic ate that the outflow s are located near the 
emission source (|Hamann et al.ll"l997r i. 

In the current work, we include all the data from our 
BAL monitoring programme. While earlier work on BAL 
variability have data on time-scales ^1 yr jBarlowl Il993l : 



Throughout this paper, all time intervals are measured in the 
rest frame of the quasar. 



ILundgren et~ai1l2007l ; IGibson et alj|2010l ; Papers 1 and 2), 
information on variability on time-scales <1— 2 months is 
very limited. Therefore, we extended our monitoring pro- 
gramme to specifically probe the shortest time-scales, down 
to ~1 week in the quasar rest-frame. As mentioned above, 
variability on short time-scales, or lack thereof, will help con- 
strain the location of these outflows, which is crucial for un- 
derstanding the physics of the flows. Furthermore, no matter 
what causes the variability, variability on the shortest time- 
scales would require extreme circumstances, such as very 
rapid continuum flux changes, very rapid crossing speeds, 
and/or small outflow structures. 

Besides examining the shortest time-scales, we also use 
the multi-epoch nature of our dataset, which spans a wide 
range of time intervals (At) from 0.02 yr (7 days) to 8.7 
yr, to investigate the incidence of C IV BAL variability as a 
function of At. We look for any drop-off in the incidence of 
variability at the shortest time-scales, and at longer, multi- 
year time-scales, we investigate whether the variability frac- 
tion continues to increase as the time-scale increases, as in 
Paper 1. 

We give an overview of the BAL monitoring programme 
in Section 12.11 and describe our analysis in Section 12.21 In 
Section 13.11 we investigate how changes in the absorption 
strength of BALs correlate with time-scale. We highlight 
specific cases of BAL variability on the shortest time-scales 
we probe in Section [3721 and we examine trends in C IV BAL 
variability with time-scale in Section 13.31 In Section 0] we 
summarize the main results of this series of papers on BAL 
variability. We then discuss, in Section 5, potential causes 
of the observed BAL variability and the implications of the 
results of this work in terms of physical properties of the 
outflowing gas. 



2 DATA AND ANALYSIS 

2.1 Observations and quasar sample 

This paper continues our analysis of the same 24 BAL 
quasars studied in Papers 1 and 2 (see table 1 in Papers 
1 and 2 for more inform a tion on the sample). This sam- 
ple was defined in IBarlowl (|l993D . which contains data from 
the Lick Observatory 3-m Shane Telescope, using the Kast 
spectrograph. We include data from this study that have a 
resolution of R = X/AX w 1300 (230 km s" 1 ), or R w 600 
(530 km s" 1 ) if there is no higher resolution data available 
for that epoch. Either resolution is sufficient for studying the 
variability in BALs since they are defined to have a width 
of at least 2000 km s" 1 . 

We have been monitoring this sample of BAL quasars 
using the MDM Observatory 2.4-m Hiltner telescope with 
the CCDS spectrograph, with a resolution of R w 1200 (250 
km s" 1 ). We also include Sloan Digital Sky Survey (SDSS) 
Data Release 6 spectra, which ha ve a resolution of R ~ 2000 
(150 km s _1 ), when available |Adelman-McCarthv et al] 
2008). We include in this work all of the observations de- 
scribed in Papers 1 and 2, and more details on the telescope 
and spectrograph setups can be found therein. 

For the current work, we supplement the data from Pa- 
pers 1 and 2 with additional observations of a subsample of 
our BAL quasars, covering rest-frame time intervals of At 
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~1 week to 1 month, during the first half of 2010. We chose 
objects with a right ascension of 08—15 hrs so that we could 
observe them from January through May. This restricts our 
sample to 17 out of 24 quasars, and out of these 17, we 
were able to observe 13 of the quasars 2 to 4 times each. 
We tended to give preference to quasars that were known to 
have variable C IV BALs in order to achieve our goal of test- 
ing for small variability time-scales in a reasonable amount 
of telescope time. Therefore, in the analysis below, we might 
be slightly biased towards higher variability fractions, espe- 
cially at the shorter time-scales. We address this bias below, 
in Section T3. 31 and find that the bias is minimal. 

Most of these observations were taken at the MDM 2.4- 
m in 2010 January (2010.04), February (2010.13), March 
(2010.20), and May (2010.35). We also have some observa- 
tions taken at the KPNO 2.1-m in 2010 February (2010.11). 
We used the GoldCam spectrograph with a 600 groove per 
mm grating in first order and a 2 arcsec slit, providing a 
spectral resolution of R « 1100 (275 km s _1 ). The observed 
wavelength coverage is ~3600 to 6100 A for all the observa- 
tions, which covers the full Siiv and C IV absorption regions 
for the quasars observed. The co-added exposure times were 
2.5—3 h per source. The typical noise level for our observa- 
tions are shown in fig. 1 of both Papers 1 and 2. 

We incorporated this new data into our existing dataset 
to create one large dataset for our sample of 24 BAL quasars, 
covering time-scales, At, of 0.02 to 8.7 yr. Table [T] summa- 
rizes the full dataset from our BAL monitoring programme, 
including the emission redshift, z em Q and the 'balnicity in- 
dex' (BI) for each object (as calculated in Paper 1). The 
fourth column gives the total number of observations, and 
the fifth column gives the range in At values for all the mea- 
surements of each object, where one measurement is a pair 
of observations. The addition of the most recent data de- 
scribed above to the existing sample gives a total of 17 BAL 
quasars with measurements that cover time-scales <0.20 yr 
(72 days). The sixth column lists the number of measure- 
ments at these short time-scales for these 17 quasars, and 
the final column gives the range in At values, in days, for 
these measurements. 



2.2 Measuring BALs and their variability 

This paper focuses on variability in C IV BALs over different 
time-scales. In order to measure variability, we first adopt 
the velocity intervals over which C IV BAL absorption occurs 
in each quasar defined in Paper 1. We used the definition of 
BI as a guide for defining these regions, i.e. they must con- 
tain contiguous absorption that reaches >10 per cent below 
the continuum across at least 2000 km s _1 . 

In order to search for variability between two epochs, 
we used a simple vertical scaling, as described in more detail 
in Papers 1 and 2, to match all the spectra to the fiducial 
Lick spectrum defined in Paper 1. We consider all available 
unabsorbed spectral regions when scaling the spectra, and 
we focus especially on the spectral regions that immediately 

2 The values of z em were obtained from the NASA /IPAC Extra- 
galactic Database (NED), which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, under con- 
tract with the National Aeronautics and Space Administration. 



border the absorption lines of interest. In some cases, a sim- 
ple scaling does not produce a good match, and there were 
disparities in the overall spectral shape between the com- 
parison spectra. For these cases, we typically fit a linear 
function, or rarely, a quadratic function, to the ratio of the 
two spectra across regions that avoid the BALs and then 
multiplied this function by the MDM or KPNO spectrum 
to match the fiducial Lick data. 

For this work, we use precisely the same criteria as de- 
scribed in Papers 1 and 2 for determining whether a quasar 
has a variable BAL. This process includes visually inspect- 
ing each pair of spectra for each quasar. For a pair of ob- 
servations to qualify as an incidence of variability, the vary- 
ing region(s) must meet two thresholds. First, the candidate 
varying region must be at least 1200 km s _1 in width. Sec- 
ond, we use the average flux and associated error within that 
velocity interval to place an error on the flux differences be- 
tween the two spectra (see Equation 1 in Paper 2). The flux 
difference in this interval must be at least 4<r to be included 
as a varying region. 

In general, we take a conservative approach and omit 
a small number of ambiguous cases that meet this thresh- 
old. Our goal here is to avoid over-counting the occurrences 
of variability. Photon statistics alone are not sufficient for 
defining real variability because flux calibrations, a poorly 
constrained continuum placement, and underlying emission- 
line variability can all add additional uncertainty not cap- 
tured by photon statistics. However, on the shortest time- 
scales, we find that the amplitudes and velo city intervals 
of va riability both tend to be small (Paper 1; iGibson et al.l 
120081 ; and, Section 3.1 below). We therefore consider two 
categories of short-term variability in our analysis, namely, 

1) "secure" cases that meet all of our criteria above, and 

2) "tentative" cases, which are likely real variations but are 
excluded by our conservative approach. We discuss several 
examples of secure and tentative detections of variability 
in Section 13.21 below. We also looked for very short-term 
changes by comparing exposures taken on different nights 
within the same observing run. The constraints are poor 
due to lower signal-to-noise in individual exposures, but we 
checked for such changes within the 2010 data and found no 
variability. 

In order to measure the absorption strength of indi- 
vidual troughs, we first adopt the power-law continuum fits 
defined in Papers 1 and 2. A power-law continuum is fit to 
one fiducial Lick spectrum per object. Since the broad emis- 
sion lines (BELs) can also vary, we only apply the power-law 
fit to the other epochs included here, and we do not fit the 
emission lines. Any measurements of absorption strength in 
this paper, as in Section[37l] are performed at velocities that 
are mostly unaffected by the BELs, so fitting the BELs for 
each epoch is not necessary. 

As in Papers 1 and 2, we first calculate the absorption 
strength, A, which is the fraction of the normalized con- 
tinuum flux removed by absorption (0 ^ A $C 1). For each 
measurement of variability, we then calculate the average 
A for each epoch within the velocity interval that varied 
and determine the change in strength, AA, between the two 
epochs. We do not divide the variable intervals into bins as 
we do in Papers 1 and 2; instead, we calculate the aver- 
age A A over the entire velocity interval that varied. In Pa- 
pers 1 and 2, we find that variability tends to occur in just 
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Table 1. Full BAL Quasar Dataset and Shortest Time-Scale Data 



Name 


z 


BI 


No. of 


At 


No. of Meas. 


At < 72 








Obs. 




<72 days 


(davs) 


0019+0107 


2.13 


2290 


7 


0.08-5.79 


2 


28-30 


0043+0048 


2.14 


4330 


5 


0.35-6.13 






0119+0310 


2.09 


6070 


3 


0.65-5.57 






0146+0142 


2.91 


5780 


1 


0.52-5.15 






0226-1024 


2.25 


7770 


2 


4.66 






0302+1705 


2.89 





3 


0.27-4.42 






0842+3431 


2.15 


4430 


13 


0.03-6.51 


9 


10-64 


0846+1540 


2.93 





8 


0.04-4.93 


4 


14-65 


0903+1734 


2.77 


10700 


7 


0.04-5.38 


2 


15-32 


0932+5006 


1.93 


7920 


13 


0.02-7.37 


9 


9-39 


0946+3009 


1.22 


5550 


9 


0.03-8.67 


7 


12-51 


0957-0535 


1.81 


2670 


8 


0.02-6.93 


7 


9-42 


1011+0906 


2.27 


6100 


10 


0.03-6.55 


5 


10-37 


1232+1325 


2.36 


11000 


3 


0.35-5.93 






1246-0542 


2.24 


4810 


7 


0.02-6.52 


3 


8-25 


1303+3048 


1.77 


1390 


9 


0.03-6.51 


7 


9-41 


1309-0536 


2.22 


4690 


7 


0.02-6.54 


3 


8-25 


1331-0108 


1.88 


10400 


7 


0.02-6.64 


1 


9 


1336+1335 


2.45 


7120 


4 


0.07-5.79 


1 


27 


1413+1143 


2.56 


6810 


7 


0.03-5.89 


1 


9 


1423+5000 


2.25 


3060 


8 


0.02-6.49 


3 


8-25 


1435+5005 


1.59 


11500 


7 


0.03-8.15 


3 


10-31 


1524+5147 


2.88 


1810 


9 


0.02-5.39 


3 


7-32 


2225-0534 


1.98 


7920 


3 


0.27-0.73 







porti ons of troughs, and sometimes in very narrow portions 
falso lGibson et al"1l200cf ). These measurements of A and A A 
give a direct measurement of the strength of the lines, and 
the change in strength, in only the portions that varied. 
Equivalent width (EW) measurements of a wide trough are 
less sensitive to strength changes in a narrow portion of the 
trough. 



3 RESULTS 

3.1 Amplitude of Variability 

In Paper 1, we directly compared the C IV BAL variabil- 
ity between a "short-term" time interval of 0.35—0.75 yr 
and a "long-term" interval of 3.8—7.7 yr, and the change 
in strength was generally slightly higher in the "long-term" 
data. To further investigate whether there is truly a corre- 
lation between the amplitude of strength changes and the 
variability time-scale, we plot AA versus At (Fig. [TJ. In or- 
der to avoid the added uncertainty of overlapping emission 
lines, we calculate AA only for absorption between —27 000 
and —8 000 km s _1 (Paper 1). We find an increasing upper 
envelope of values of A A with increas i ng A t. This is consis- 
tent with the findings of lGibson et al. (2008), who show that 
the envelope of values of AEW expands as At increases. We 
only plot variable cases in Fig. 1, as well as measurements 
of the tentative cases of variability (green triangles; Section 
12. 2|) . On the shortest time-scales that we are probing in this 
paper, with At < 0.20, the AA values are all below ~0.1. 
On multi-year time-scales, the values of AA reach as high 
as ~0.46. 



0.5 r 




o.o L_i , , , , , - 

0.1 1.0 10.0 

At (yeors) 



Figure 1. The change in absorption strength, AA, for variable 
portions of absorption troughs versus the time interval between 
observations, At. The green triangles at shorter time-scales rep- 
resent measurements of AA for the tentative cases of variability. 

3.2 Examples of Variability on the Shortest 
Time-Scales 

As mentioned in Section [21 we obtained new data specifi- 
cally to augment the temporal sampling of our dataset at 
time-scales <0.20 yr (<~72 days) in the quasar rest-frame. 
We specifically look at this range of time-scales because it 
corresponds to the four shortest time-scale bins in Fig. 9 
below, where we examine trends with the BAL variability 
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Figure 3. Spectra of just the Civ BAL in 1246-0542 for the 
same 3 epochs shown in Fig. [2] The shaded regions mark the 
intervals of variability. As in Fig. [2] the black dashed curve shows 
the power-law fit to the continuum, and the formal la errors for 
each spectrum are shown near the bottom. 

time-scale. Out of the 17 quasars for which we have data 
on these short time-scales, only 2 quasars exhibited C IV 
BAL variability (12 per cent). If we include the tentative 
cases of variability, then 5 out of these 17 quasars varied 
(29 per cent). For comparison, even with the tentative cases 
included, this fraction is still lower than the incidence of 
variability in our "short-term" sub-sample from Paper 1, 
where 39 per cent of the quasars varied over time intervals 
of 0.35-0.75 yr. 

In this section, we first highlight the only two quasars 
with secure detections of C IV BAL variability on the short- 
est time-scales in our data set (<0.20 yr). We then describe 
examples of cases that were classified as tentative detections 
of variability on these short time-scales. We emphasize that 
all of the cases of variability, and tentative variability, dis- 
cussed here meet the threshold defined in Paper 1 and used 
throughout this paper, i.e. the variability must occur over 
a region at least 1200 km s" 1 wide and the flux difference 
between the two spectra must be at least 4a (see Section 
I2.2|l . We present these individual cases in order of increas- 
ing time-scale. 

3.2.1 1246—0542: Secure variability over 8 Days 

The quasar 1246—0542 is one of just two quasars in our sam- 
ple for which we have a secure detection of variability over 
a time-scale <0.20 yr. Our dataset contains measurements 
down to At of 0.02 yr (~7— 8 days), and we found variabil- 
ity down to this time-scale in 1246—0542. We present this 
measurement in Fig. [5J showing the full spectrum for the 
2010.13 and 2010.20 observations, and in Fig. [3] showing 
just the C IV BAL. As indicated by shaded bars in Figs [2] 
and [3] there are two separate velocity intervals of variability 
between the 2010.13 (black curve) and 2010.20 (red curve) 
observations. The two velocity intervals are centered at —24 
050 km s _1 and -17 200 km s _1 , with widths of 1300 km 
s _1 and 2600 km s _1 , respectively. The flux differences in 
these two intervals are 8.0a and 13a. 

The rightmost shaded bars in Fig. [2] mark the variable 



intervals in C IV, and they correspond to the shaded inter- 
vals in Fig. [3] We shifted these intervals to the corresponding 
velocities in Nv, marked by the leftmost shaded bars, and 
Si IV. The doublet separations for N v and Si IV are 960 and 
1930 km s~ , respectively, which are both wider than the 
doublet separation for C IV (500 km s _1 ). When shifting the 
shaded intervals to Nv and Siiv, we therefore widened the 
intervals based on the greater doublet separations in these 
lines. The identification of N v absorption, and any potential 
variability, is complicated by Lya forest absorption, which 
appears throughout the entire spectrum blueward of the Lya 
A1216 emission line. There is also a slight mismatch between 
the two spectra near the limit of the wavelength coverage. 
However, we use the velocities of C IV absorption to identify 
where the N V absorption should be located, and the veloc- 
ities of apparent variability in Nv, in the redmost velocity 
interval especially, match the velocities at which C iv ap- 
pears to vary. Corresponding variability in N v supports the 
detection of variability in C iv. We also note that despite 
the slight mismatch at the limit of the wavelength cover- 
age, the two spectra match well along the remainder of the 
spectral coverage and, in particular, in the spectral regions 
that immediately border the C iv absorption. We do not find 
any corresponding changes in the Siiv BAL. In Paper 2, we 
found just one tentative case where C iv variability was not 
accompanied by Si iv variability at the same velocity. 

Figs and also show that the BAL in 1246-0542 
again varied between 2010.20 and 2010.35 (At ~ 0.05 yr, 
or 17 days), returning to its previous state from 25 days 
earlier. There is no difference in the BAL strength between 
the earlier 2010.13 spectrum (black curve) and the 2010.35 
spectrum (blue curve). We also have measurements of this 
quasar at other epochs with variability in the same velocity 
intervals within which we detect variability in Figs [2] and |3j 
In particular, in Papers 1 and 2, we find variability, with a 
flux difference of 6a, between the two "long-term" epochs, 
1992.19 and 2008.35, at the same velocities as the bluemost 
interval here (see fig. 1 in Papers 1 and 2). These "long- 
term" epochs are two different observations that were taken 
at earlier epochs than the observations shown here in Figs 
[3] and [3] This gives further confirmation that the detection 
of variability over 8 days between 2010.13 and 2010.20 is 
secure and not a result of a poor flux calibration in one of 
the spectra or a poor match between the unabsorbed regions 
in the two spectra. 

3.2.2 0842+3431: Secure variability over 10 Days 

The quasar 0842+3431 is the only quasar besides 1246-0542 
with a secure detection of variability at At < 0.20 yr. We 
detect variability over 0.03 yr (10 days) between the 2010.11 
and 2010.20 observations of 0842+3431. These spectra are 
shown in the top panel of both Figs|4]and[5] and the variable 
interval, marked by a shaded bar, is centered at —17 800 km 
s _1 and is 1800 km s" 1 wide. The flux difference between 
the two epochs is 8.5a. 

This variable interval is on the blue side of the trough 
at a high enough velocity that this detection of variability is 
not affected by any variability in the Civ BEL. The spectra 
also match very well on either side of the shaded interval. 
Furthermore, we mark the corresponding velocities in Siiv 
with a shaded bar, and there is potential variability there. 
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Figure 2. The full spectra of 1246—0542, showing 3 epochs covering a At of 0.02—0.07 yr. The emission lines are identified by vertical 
bars. The rightmost horizontal bar marks the Civ BAL absorption, and the shaded regions below it mark the intervals of variability. 
The other horizontal bars and shaded regions mark the corresponding velocities in Nv and Si IV. We use binomial smoothing to improve 
the presentation of these and all of the other spectra displayed in this paper. The black dashed curve shows the power-law fit to the 
continuum, and the formal lcr errors for each spectrum are shown near the bottom. 



The middle panel of Fig. U shows two Lick spectra, 
1990.90 and 1991.11, which are separated by 24 days, and 
the bottom panel shows the 1991.31 and 1991.86 spectra, 
which are separated by 64 days. In both of these measure- 
ments, there is variability in Siiv at the same velocities 
as the variability in Civ. Between 1990.90 and 1991.11, 
the variability is stronger towards the blue and red ends 
of the interval, but much weaker in the middle of the in- 
terval. The variability in Siiv, however, is most significant 
in the middle of the shaded interval. All four spectra from 
1990.90 to 1991.86 are shown in the bottom panel of Fig. [5] 
(a full discussion of the multi-epoch nature of the variabil- 
ity in 0842+3431 can be found in Paper 2). The clear low- 
amplitude variability on 24 and 64-day time-scales in the 
Lick data support our determination of similar variability 
(albeit at different velocities) over 10 days between 2010.11 
and 2010.20. 

3.2.3 1011+0906: Tentative variability over 17 days 

In this section, we describe a measurement that we classified 
as a tentative detection of variability. The full spectra for the 
2010.20 and 2010.35 observations of 1011+0906, which are 
separated by a At of 0.05 yr (17 days), are plotted in Fig. 
[6] showing two potential regions of variability with shaded 
bars. These intervals are centered at —26 250 and —20 400 
km s _1 and are 1900 and 1400 km s _1 wide, respectively. 
The significance of the variability is 9.1a in the wider inter- 
val and 6.9a in the narrower interval. These values are well 



above the threshold of 4a, and there is potential variability 
in Siiv at the same velocities as in Civ. However, there ap- 
pear to be mismatches throughout the spectrum, blueward 
of the Siiv emission line. These mismatches are similar in 
amplitude to the amplitude of the strength changes in the 
regions of potential C IV variability identified by the shaded 
regions in Fig. [HJ We considered the possibility that the slope 
of one of the spectra needs to be adjusted to better match 
the other. We examined the other data collected on the same 
nights as these two spectra, but found no evidence for any 
issues with the flux calibration at either epoch. However, 
we maintain our conservative approach here and keep this 
measurement in the "tentative" category. 

3.2.4 0932+5006: Tentative variability over 31 days 

We present here another example of a tentative detection 
of variability. Fig. [7] shows the C IV and Si IV BALs for the 
2008.03 and 2008.28 observations of 0932+5006, which are 
separated by 0.09 yr (31 days). The shaded region identifies 
the candidate variable interval, with a width of 1200 km s _1 , 
which is just at the threshold defined in Paper 1. The flux 
difference within this interval is 8.3a, and there is also pos- 
sible variability at these velocities in Si IV. As in the case of 
1011+0906, however, there are mismatches in other regions 
of the spectra, in particular towards the blue end of the spec- 
tra in the bottom panel of Fig. [7] Given the small amplitude 
of variability in an interval that is just wide enough to meet 
our threshold, combined with mismatches in other regions 
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Figure 4. Spectra of 0842+3431. Each panel shows a measurement of variability over a time-scale of <0.20 yr, with the top panel 
showing two epochs separated by just 0.03 yr (10 days). The rightmost shaded regions mark the intervals of Civ variability for each 
measurement, and the leftmost shaded regions mark the corresponding velocities in Si IV. 



of the spectra, we classify this measurement as a tentative 
case of variability. 



3.2.5 0903+1734: Secure variability over 72 days 

The variability between 2006.31 and 2007.05 in 0903+1734 
is one of the most prominent cases of variability on time- 
scales less than those analyzed in Paper 1, i.e., <0.35 yr. 
These observations are separated by 0.20 yr (72 days), which 
is just on the upper threshold of the range of time-scales 
included in the fourth point of Fig. [9] below. Fig. [8] shows 
these two spectra with the variable intervals in C IV, and the 



corresponding velocities in Si IV, shaded. The variability here 
covers a wide range in velocities, unlike the other examples 
at shorter time-scales in this section. Paper 2 contains a full 
discussion on the multi-epoch behavior of this quasar. 



3.3 Time-scales for Variability 

In this section, we examine the relationship of C IV BAL 
variability with time-scale across the full measured range 
from 0.02 to 8.7 yr. To do this, we compare the BALs in 
each pair of observations at all velocities in every quasar 
and then count the occurrences of C IV BAL variability, us- 
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Figure 6. The full spectra of 1011+0906, showing 2 epochs covering a At of 0.05 yr. We classify this measurement as a tentative case of 
variability. The rightmost shaded regions mark the candidate intervals of C IV variability, and the leftmost pair of shaded regions mark 
the corresponding velocities in Si IV. 
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Figure 8. The full spectra of 0903+1734, showing 2 epochs with a At of 0.20 yr. The four rightmost shaded regions mark the intervals 
of Civ variability, and the leftmost shaded regions mark the corresponding velocities in Si IV. Strong narrow lines located at ^1425 A 
and blueward of 1300 A are due to ions at much lower redshifts, unrelated to the quasar. 
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Figure 5. Spectra of the Civ BAL in 0842+3431. The top panel 
shows the same two epochs from the top panel of Fig. [4] and the 
bottom panel shows all four Lick epochs from the bottom two 
panels of Fig. [4] The shaded regions mark the extent of the C IV 
variability. 




-3x10" -2x10" -1x10" 

Velocity (km/s) 

Figure 7. Spectra of the C IV (top panel) and Si IV (bottom panel) 
BALs in 0932+5006 for two epochs separated by 0.09 yr (31 days). 
This is another tentative measurement of variability. 



ing our definition of BAL variability first defined in Paper 1 
(see Section f2 . 2 [) . We then calculate a probability by divid- 
ing the number of occurrences of variability by the number 
of measurements, where a pair of observations is one mea- 
surement, in logarithmic bins of At. We plot this measured 
probability of detecting Civ BAL variability versus At in 
Fig. [9] The la error bars are based on counting statistics 
for the number of occurrences of variability a nd the number 
of measurements in each bin |Cameronll201ll ). 

Fig. [9] confirms our results from Paper 1 that the inci- 
dence of C IV BAL variability is higher when At is larger. 
The probability of C IV variability approaches unity when At 
^ 10 yrs. This decreases to ~0.6 for At ~1 yr and ~0.05— 0.1 
for At ~0.08 yr (~1 month). Below At ~ 0.08 yr, the mea- 
sured probability of variability begins to increase slightly as 
At becomes shorter. For the four shortest time-scale bins, 
where At <0.2 yr, we also calculate this measured proba- 
bility of variability with the tentative cases of variability, as 
described in Sections 12.21 and 13.21 included. The resulting 
values are plotted as green triangles in Fig. [5] This tends 
to flatten out the trend in the shortest time-scale bins. The 
green and black dotted curves show the least-squares fit to 
the points with and without the tentative cases of variabil- 
ity included. Overall, inclusion of the tentative cases of vari- 
ability does not significantly change the slope of the points. 
We note that these are probabilities for detecting variabil- 
ity between two measurements separated by At, and not for 
detecting variability at any time in that At time-frame. In 



other words, if a quasar varied then returned to its initial 
state all within a certain At, then a measurement at that 
value of At would not count as variable in this plot. 

The values plotted in Fig.[!TJrepresent probabilities only 
if each measurement is an independent event. There are two 
primary biases that affect this plot. First, some quasars con- 
tribute more than others (see Table 1), and those quasars 
that were observed most frequently throughout our monitor- 
ing programme were typically ones that were known to be 
variable (Section l2.ip . Second, our observations are clustered 
in time, with most taken at Lick from 1988 to 1992 and at 
MDM from 2007 to 2010. The middle and bottom panels of 
Fig. [9] display the number of measurements and the number 
of quasars, respectively, that contribute to each At bin (each 
quasar can contribute multiple times to a single bin). The 
middle panel shows that the majority of the measurements 
are clustered at shorter (<1 yr) and longer (>~4 yr) time- 
scales. The average number of measurements per quasar in 
each bin ranges from TO (in the shortest time-scale bin) to 
8.5 (in the penultimate bin). 

Since those quasars contributing the most measure- 
ments to Fig. [9] are typically those that were known to be 
variable, the measured probabilities of variability in Fig. [9] 
may be biased to higher values. Those quasars with the most 
epochs of data will be counted much more often in Fig. [9] 
For example, we have 13 epochs of data for 0842+3431 and 
0932+5006, and up to just 10 epochs of data for each of 
the other quasars in the sample. Thirteen epochs for one 
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Figure 9. The top panel displays the fraction of measurements 
in which Civ BAL variability was detected at >4<r over a velocity 
interval of at least 1200 km s _1 versus the time interval between 
the two observations. The green triangles show the variability 
fractions with tentative cases of variability included. The middle 
and bottom panels show the number of measurements and the 
number of quasars, respectively, that contribute to each bin. 



quasar provides 78 pairs of observations, so 0842+3431 and 
0932+5006 together contribute 27 per cent of the measure- 
ments in Fig. [9] To address this bias, we therefore show 
in Fig. [TO] a version of this plot which omits these two 
quasars with the most data. We plot a least-squares fit to 
the data points in Fig \W\ (dashed curve) and overplot the 
least-squares fit from Fig. (dotted curve) to compare their 
slopes. This shows that the most frequently observed quasars 
are not significantly affecting the slope of the points, but 
they do shift the line to slightly higher probabilities. We do 
not include the tentative cases of variability in this figure, 
or in Figs 1111 through 1131 below. 

Even though the slope does not change significantly be- 
tween Figs 191 and 1 101 there is one point at At of 2.4 yr that 
does decrease significantly in Fig. [10] In Fig. [9] there are only 
6 quasars that contribute measurements to this bin because 
a spectrum from SDSS is necessary for having a measure- 
ment at these intermediate time-scales. We only have SDSS 
spectra for 8 out of the 24 quasars in our sample. By re- 
moving 0842+3431 and 0932+5006, we are left with just 4 
quasars in this bin centered at 2.4 yr, and more than half 
of the measurements of variability in this bin in Fig. [5] are 
measurements of 0842+3431 and 0932+5006. 

This relates to the second bias mentioned above that 
instead of randomly sampling the quasar light curves, our 
measurements are clustered at shorter (<1 yr) and longer 
(>~4 yr) time-scales. For example, in cases where we have 
just Lick and MDM/KPNO data, any variability that we 
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Figure 10. Same as Fig. [9] but with the two most-observed 
quasars, 0842+3431 and 0932+5006, removed. The dashed line 
in the top panel is a fit to the data points, and the fit from Fig. 
[9] is overplotted as a dotted line. 




At (years) 

Figure 11. A version of Fig.[9]where each quasar only contributes 
once to each At bin, and then we iterate 1000 times to get a range 
of least-square fits, shown here by the shaded region. The dotted 
line is the least-squares fit from Fig. [9] 
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detect between the last Lick observation and the first MDM 
observation for that quasar could have occurred within any 
time interval during those ~4 years. In this scenario, if there 
were, for example, 2 Lick observations and 5 MDM obser- 
vations for an object, and the only variability that occurred 
in this quasar happened between the last Lick and the first 
MDM observations, then this quasar would contribute 10 
measurements of variability to the last two At bins. We 
therefore may be biased towards greater variability fractions 
at longer time intervals over intermediate and short time in- 
tervals. 

To explore any bias that our uneven sampling may in- 
troduce into Fig. [9l we created a version of this plot where 
each quasar can only contribute one measurement to each 
At bin. If a quasar has multiple measurements in a par- 
ticular bin, then one of those measurements is chosen at 
random to be included in this plot. We then recalculate this 
plot 1000 times, and for each iteration, we randomly choose 
which measurement per object is included in each bin and 
then create a least-squares fit to the resulting probability 
values. In Fig. 111! we show the range of least-squares fits 
we obtain. By iterating many times, we are attempting to 
uncover a range of slopes that includes what the slope would 
be in the ideal case of even sampling across the entire range 
of time-scales covered here. The range in slopes that we ob- 
tain overlaps the slope of the points in Fig. [5] shown here 
again by the dotted curve, but the slopes tend to be shal- 
lower and/or shifted to slightly lower probabilities across 
the full range in At. This indicates that we do indeed have 
a slight bias towards greater variability fractions at longer 
time-scales in Fig. [9] At the shortest time-scales, the curve 
from Fig.[5]is roughly in the middle of the range of curves in 
Fig. 1111 We mention above in Section \2. II that we tended to 
monitor those quasars known to vary, which might introduce 
a bias towards greater variability fractions at the shortest 
time-scales. However, it is not clear from Fig. [11] whether 
such a bias exists, but if it does, it is certainly weaker than 
any bias at longer time-scales. 

Finally, we present Fig. 1121 which shows the measured 
probability of C IV BAL variability versus time-scale, consid- 
ering only the velocities that contribute to the measurement 
of BI, i.e., -25 000 to -3 000 km s" 1 . Therefore, we only in- 
clude absorption i n this plot that are BAL s according to the 
strict definition of IWevmann et all (|l99ll ). so the 2 quasars 
in our sample with BI=0, 0302+1705 and 0846+1540, are 
not included at all here. Restricting to this velocity range 
also removes uncertainty in the variability detection due to 
the presence of a potentially variable broad emission line 
(BEL). 

The dotted line in Fig. I12l represents the slope from Fig. 
[5] and the dashed line is a fit to the points in the current fig- 
ure. These lines show that the restriction in velocity range 
changes the slope of the points; the slope becomes shallower. 
In Paper 1, we found that the incidence of variability in- 
creases at higher outflow speeds. Here, we are removing the 
highest-velocity measurements, so that will shift the points 
to lower values. If the incidence of variability at velocities 
<— 25 000 km s _1 is the same across the entire range of 
time-scales, then the points at longer time-scales will be af- 
fected more than the points at smaller time-scales, where 
the overall incidence of variability is smaller. For example, 
if excluding the highest velocities lowers the recorded inci- 
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Figure 12. Same as Fig. [9] but here we only consider absorption 
that lies between —25 000 and —3 000 km 8 , thus satisfying all 
of the requirements of the balnicity index. The dotted line is the 
fit to the data points in Fig. [9] and the dashed line is a fit to the 
data points in the current figure. 



dence of variability by 20 per cent at all time-scales, then 
the points at short-timescales, where the measured probabil- 
ity of CIV BAL variability is small, will only decrease very 
slightly, wherease the points at longer time-scales, where the 
probabilities are greater, would decrease by nearly 0.2. Thus, 
the slope in Fig.[l2]would naturally be shallower than in Fig. 

El 

We also investigate the cumulative probability of Civ 
BAL variability. In Fig. 1131 each point represents the frac- 
tion of BAL quasars with variable absorption at any At 
value up to that point. We calculate this distribution using 
the full measurement sample in Fig. [9] except that we omit 
2225—0534 because we only have data over time-scales of 
At <1 yr for this object. By the largest At values, all the 
quasars in our sample, except for 2225—0534, are contribut- 
ing. As in Fig. 1111 we only allow each quasar to contribute 
once to each bin, and we then iterate 1000 times to obtain 
a range of variability fractions for each bin. This range is 
represented by the shaded region in Fig. 1131 where the black 
histograms indicate the minimum and maximum variability 
fractions at each time-scale. If we included all measurements 
of each quasar in each bin, then as long as a quasar varied at 
least once in a particular bin, it would be included as vari- 
able in that bin and in all of the higher time-scale bins. In 
this case, the resulting curve would be the upper histogram 
in Fig. [H 

There is no range of values for the first four bins of 
Fig. [T3] because of the limited number of measurements per 
quasar in this time-scale regime. Below ~0.2 yr, the fraction 
of quasars exhibiting Civ BAL variability is ~13 per cent. 
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Figure 13. A cumulative distribution of the fraction of quasars 
with C IV variability with time-scale. This distribution is drawn 
from the full data set shown in Fig. [9] excluding 2225—0534. As 
m Fig. HU each quasar only contributes once to each bin, and 
we iterate 1000 times to obtain a range of variability fractions at 
each time-scale. Because of the small number of measurements at 
the shortest time-scales, there is no range of values in the first 
four bins. 



At time-scales up to ~1 yr, this fraction increases to 45 ± 
23 per cent, and over multi-year time-scales, the variability 
fraction reaches 83 ± 9 per cent. 



4 SUMMARY OF RESULTS 

This work is the third paper in a series on BAL variability 
in a sample of 24 BAL quasars, with observations covering 
a wide range of rest-frame time-scales from 0.02 to 8.7 yr. 
Paper 1 describes general trends in the C IV BAL variability 
data and finds, in particular, that variability occurs more of- 
ten at higher outflow velocities and in shallower troughs (see 
alsolLundgren et al.1120071 ; iGibson et ai1l2008l: iFiliz Ak et al] 
l2012n . In both Paper 1 and the multi-epoch analysis of Paper 
2, we note that variability typically occurs in just po rtions 
of troughs (e.g., Figs [2j— [5l see also lGibson et alJl200Sh . In 
rare cases, BAL features appear, dis appear, or change to 
or from narrower mini-BAL features l|Hamann et al.l 120081 ; 
Leighlv et alll2009l; iKrongold et , al.l l2010l : Papers 1 and 2; 



Vivek et al.ll2012i ; lFiliz Ak et al.ll20L?T ) 



Paper 2 also directly compares variability in C IV ab- 
sorption to variability in Siiv and determines that Siiv 
BALs are more likely to vary than Civ BALs. For exam- 
ple, at flow speeds >— 20 000 km s _1 in the "long-term" 
sample defined in Paper 1, 47 per cent of quasars exhib- 
ited Siiv variability while only 32 per cent exhibited ClV 
variability. Furthermore, approximately half of the variable 
Si IV regions did not have corresponding C IV variability at 
the same velocities, while in just one poorly measured case 



were changes in C IV not matched by variability in Si IV. The 
greater variability in Siiv is likely due to the tendency for 
weaker lines to vary more (Paper 1) , combined with the fact 
that the Si IV absorption is typically weaker than C IV (Pa- 
per 2). Furthermore, when Civ and Siiv both varied at the 
same velocity, the changes always occured in the same sense 
(both lines either became stronger or both weakened). Sim- 
ilarly, BAL changes at different velocities within the same 
line almost always occurred in the same sense. 

In the present paper, we examine C IV BAL variabil- 
ity behaviours as a function of time-scale, with a particular 
emphasis on the shortest times. Over 50 per cent of mea- 
surements separated by 1 yr in the quasar rest-frame varied, 
and on multi-year time-scales, the probability of detecting 
variability between two observations of the same quasar ap- 
proaches unity (Fig. [9} . We therefore detect a strong trend 
towards greater variability fractions over longer time-scales, 
which is consistent with our findings in Paper 1. 

Not only does the incidence of variability decrease at 
shorter time-scales, but the amplitude of the absorption 
strength changes (A A) also decreases (Fig. [TJ . Over multi- 
year time-scales, we find a maximum AA =~0.46. However, 
all of the variations across time-scales of At < 0.27 yr have 
amplitude changes of A A < 0.2, and the shortest times we 
measure, At < 0.2 yr, have AA < 0.1. This is generally 
consistent with Gibson et all (2008), who find a decrease in 
AEW values at shorter time-scales. 

We detect variability down to the lowest At values that 
we probe (At ~ 0.02 yr, or 8-10 days; Section The 
amplitude changes on this short time-scale are small (AA 
< 0.1), and they appear in just small portions of the BAL 
troughs (in intervals of Av from 1300 to 2600 km s _1 ; Figs 
[2HS]in Sections [3T2~TH3~2~2l . 

We considered the possibility that the trend we find for 
greater variability fractions on longer time-scales may be bi- 
ased due to uneven sampling of the quasar light curves and 
a tendency to monitor quasars that already varied. As dis- 
cussed in Section [33] most of our observations are clustered 
at shorter (<1 yr) and longer (>~4 yr) time-scales. If we 
tend to monitor quasars already known to be variable, then 
those quasars will contribute many more measurements at 
At >~4 yr than those that did not vary. However, Figs 1101 
and 1111 which remove the main sampling biases, still show 
strong trends for greater variability on longer time-scales. 
This trend is also readily apparent in Fig. 1131 where we plot 
the fraction of quasars instead of the fraction of measure- 
ments. The fraction of quasars that exhibited variability at 
any At $;~0.2 yr is ~13 per cent. This fraction increases to 
45 ± 23 per cent by At~l yr and then to 83 ± 9 per cent 
over multi-year time-scales. This is greater than the value of 
65 per cent we found for the two-epoch analysis on similar 
time-scales for this s ame sample of quasa rs in Paper 1. Al- 
though, we note that IGibson et al.l (|2008r ) reported that 92 
per cent (12/13) of their quasars exhibited Civ BAL vari- 
ability in their two-epoch multi-year analysis. 



5 DISCUSSION 

5.1 Causes of Variability 

The time-scales of quasar BAL variability are important for 
constraining the location of the outflowing gas, but the spe- 
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cific constraints are dependent on the cause (s) of the vari- 
ability and the structure of the outflows. In Paper 2, we 
discuss two scenarios that could produce the observed BAL 
variability: f) changes in the far-UV continuum flux caus- 
ing global changes in the ionization of the outflowing gas, 
and 2) outflow clouds moving across our lines-of-sight to 
the quasar continuum source. We found evidence to support 
both scenarios, but for most individual cases, the cause of 
the variability is ambiguous. 

In quasars that have variable absorption at more 
than one velocity, the variabilities tend to be coordinated 
(for example, 1 246-0 542, Fig. and 0903+1734, Fig. 
lHamann et al.l |201lf ) found coordinated line variations in 
multiple narrow absorption line (NAL) systems and argue 
that this supports the scenario of a global change in ioniza- 
tion. Changes in the ionizing flux incident on the entire out- 
flow should cause global changes in the ionization of the flow. 
A change in covering fraction due to moving clouds is less 
likely when absorption regions at different velocities vary in 
concert because this would require coordinated movements 
among outflow structures at different outflow velocities and 
radii. 

On the other hand, variability in narrow por- 
tions of BAL troughs fits more naturally in the sce- 
nario of crossing clouds. This latter scenario is fa- 
vored by prev i ous work on BAL varia b ility, including 



iLundpen et al.l (120071). I Gibson et all (120081), Hamann et all 
j20Qg|). iKrongold et all (l2010h. lHall et all l|201lf) . and 



IVivek et all (|2012l ). lHamann et all l^POSl ) and Capellupo et 



al, in preparation, discuss individual cases where the C IV 
BAL is saturated. Capellupo et al, in preparation, report on 
the detection of variable Pv Af If 8, If 28 absorption at the 
same velocities as the C IV and Si IV absorption and variabil- 
ities in the same spectrum. The presence of a low-a bundance 
line l i ke Pv indicates that C iv is very saturated (|Hamannl 
1 19981 ; lHamann et al. I l2002l ; Capellupo et al, in preparation; 
and references therein). The depths of saturated lines are 
governed by the covering fractions in those lines, so a change 
in covering fraction would change the depth of the line. A 
change in ionization would produce no detectable change in 
a highly saturated absorption trough. Therefore, variabil- 
ity in saturated absorption li nes strongly favo rs the crossing 
cloud scenario. Furthermore. iMoe et alj (|2009h detect a vari- 
able BAL component in a quasar spectrum, and they con- 
clude that the variability was most likely caused by trans- 
verse motion of the outflow because they do not detect a 
significant change in the continuum flux of the quasar. How- 
ever, variability on very short time-scales might pose a prob- 
lem for the crossing cloud scenario, as we discuss below (Sec- 
tion [5X2]). 

Other possible causes for the BAL variability are more 
problematic. For example, a change in the size of the contin- 
uum source would cause a change in the covering fraction of 
saturated absorption lines. The main problem with this sce- 
nario, though, is that it should change the covering fraction 
of an entire trough and not just small portions thereof. Most 
cases of BAL varia bility occur in just portions of troughs 
jGibson et alj|2008l ; Papers I and 2). 

Another possibility is instabiliti es within th e flow s 
themselves, as in the simulations of iProga et all {2000), 
causing the observed variability. The time-scale for such 
instabilities significantly affecting the outflow structure is 



probably similar to the outflow dynamical time. If a flow 
has a radial velocity of 15 000 km s _1 and is located nomi- 
nally around the broad-emission line region (~0.I5 pc), the 
dynamical time is roughly ~10 years. This is much larger 
than many of the measurements of variability time-scales 
in this work and cannot explain the very short time-scale 
changes we observe (as in Section 13. 2[) . Furthermore, the 
time-scale for clouds dissipating without sheering or exter- 
nal forces is given by the sound crossing time. For a nominal 
cloud radius of 0.1 pc and temperature of 10 4 K, the sound 
crossing time is ~8000 yr. This is much greater than the 
variability time-scales reported in this work, so dissipation 
of the outflow clouds is most likely not causing the observed 
variability. 

As in Paper 2, we conclude here that, in general, the 
cause of variability is either changing ionization, transverse 
cloud movements, or some complex amalgam of the two. We 
therefore focus on just these two scenarios below. 



5.2 Implications of Variability 

5.2.1 Changing ionization 

Significant ionization changes in the BAL gas require at least 
a recombination time, nominally t rec ~ l/n e a r , if the gas 
maintains ionization equilibrium, where n e is the electron 
density and a r is the recombination rate coefficient. For the 
cases of 1246—0542 and 0842+3431, we can use the vari- 
ability time-scale of 8—10 days (Sections 13. 2. II and 13.2.21) to 
set a lower limit on the density of n e >^4 — 5 x 10 5 cm" 3 . 
Based on the photoionization calculations of lHamann et all 
( 1995) and this constraint on the density from the recombi- 
nation time, the maximum distance is <~200 pc, assuming 
the ionization is at least as high as that needed for a maxi- 
mum Civ/C ratio. However, the actual recombination time 
can be shorter, resulting in lower derived densit ies and larger 
distances, if the gas is mo re highly ionized |Krolik fc Krissl 
1 19951 ; lHamann et ; al.|[l997b . 

Another important implication of ionization changes 
is that they require significant variations in the quasars' 
incident (ionizing) flux. This seems unlikely for at least 
the shortest time-scales we measure because our sample 
consists of luminous quasars. More luminous quasars have 
a smaller amplitude of continuum variability than fainter 
AGN, and the amplitude of continuu m variability tends 
to de c rease on short e r tim e-scales (e.g. IVanden Berk et all 
|2004| ). iMisawa et all (|2007l ') detected variability in a mini- 
BAL over a time-scale of 16 days, and they also state that 
variability on such a short time-scale is much faster than 
any expected changes in the continuum emission of a lumi- 
nous quasar. They propose that a screen of gas, with varying 
optical depth, is located between the continuum source and 
the absorbing gas. This screen could be the ionized X-ray- 
shielding gas in the outflow models of iMurrav et all l|l995h . 
and its proposed location is just above the accretion disk. 
The screen is in co-rotation with the disk, and if the screen is 
clumpy, then as it rotates, the intensity of the ionizing con- 
tinuum that is transmitted through the screen varies. Since 
the screen is closer to the continuum source, it would be ro- 
tating at a faster velocity than the absorber and could then 
cause changes in the ionization of the absorbing gas more 
quickly than continuum variations alone. Several models in- 
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Figure 15. Transverse velocity and distance of the outflowing 
gas versus rest-frame time-scale for the two models illustrated in 
Fig. 1141 and two typical values of AA (Fig. [TJ. 



Figure 14. A simple schematic of the two models used to calcu- 
late crossing speeds of an outflow in the changing covering frac- 
tion scenario: (a) the "crossing disk" and (b) the "knife edge" 
model. The darker regions represent the outflow moving across 
the continuum source (light gray region) along our line-of-sight 
to the quasar. The arrows represent the distance the outflow must 
travel to cover the same fraction of the continuum source in both 
models a and b. 

deed predict variability on very short time-sca l es related 
to th i s X-ray shielding ga s l|Schurch et al.l 120091 : ISim et all 
l2010D . |Schurch" et al.ll2009l determine that the total column 
density in the gas could change by a factor of 2 in ~9 days. 
(However, see lHamann et al.N201ll and in preparation for 
complications related to this scenario.) Another possibility 
for producing ionization changes on very short time-scales 
is small hot spots on an inhomogeneous accretion disk (e.g. 
iDexter fc Ago] 1201 ll ) which could appear and disappear on 
short time-scales. However, all of these scenarios predict 
global changes across BAL troughs, which we typically do 
not see (Section l5.ip . 

5.2.2 Crossing clouds 

We now consider a simplified scenario where the absorbers 
have a constant ionization and column density and the vari- 
ability is due to a single outflow component moving across 
the continuum source. We can use the time-scale of the ob- 
served variability to estimate the crossing speed given a ge- 
ometry for the emission and absorption regions. We esti- 
mate a characteristic diameter for the continuum source at 
1500 A using the observed fluxes at this wavelength from 
iBarlowl i|l993l ) and a standard bolometric correction factor, 
L w 4.4AL A (1500A) (|Hamann et al.l 120111 ). in a cosmology 
with H = 71 km s" 1 Mpc, fi M = 0.3, fi A = 0.7. This yields 
bolometric luminosities of ~ 2 x 10 46 to 3 x 10 47 ergs s _1 
across our sample. Based on the average bolometric lumi- 
nosity, ~ 7 x 10 46 ergs s _1 , a characteristic diameter for the 
continuum region at 1500 A is -D1500 ~ 0.008 pc and for the 
C IV BEL region is Dciv ~ 0-3 pc, assuming L = l/3L e rfj 
and Mbh = 1.4 x 1 9 M^ (|Peterson et al.ll2004l ; lBentz et al.l 
2007; Gaskcll 2008 ; Hamann & Simon, in preparation). 



Using the characteristic diameter for the continuum re- 
gion, we calculate the crossing speed of the outflow compo- 
nent for two models, as shown in Fig. 1141 In the first model, 
a circular disk crosses a larger circular continuum source 
along a path through the center of the background circle (the 
"crossing disk" model). In the second model, the continuum 
source is square, and the moving absorber is larger, with a 
straight edge that moves across the background square (the 
"knife edge" model). To calculate the crossing speed of the 
absorber, we take the change in strength of the line in the 
varying region, AA, as the fraction of the continuum source 
that the absorber crosses within the elapsed time between 
observations, At. In the "crossing disks" model, the distance 
traveled by the absorber is then yAA D1550, and the cross- 
ing speed is this distance traveled divided by At. This model 
gives the maximum transverse velocity. In the "knife edge" 
model, the velocity is smaller than for the "crossing disk" 
model by a factor of %/ AA, and it gives the minimum trans- 
verse velocity. These two models represent opposite extremes 
which should encompass the possibilities for real absorbers. 
We plot the transverse velocity as a function of time-scale on 
the right ordinate in Fig. [TS] for a change in BAL strength, 
AA, of 0.05 and 0.2. A AA of 0.05 is a typical value for data 
on time-scales <0.1 yr, and values between AA of 0.05 and 
0.2 are typical on time-scales of 0.1—3 yr (Fig. [I}. 

Next we tie the crossing speeds to a physical location by 
assuming they are roughly equal to the Keplerian rotation 
speed around an SMBH with a mass of Mbh ~ 1.4xl0 9 Mq. 
This is a reasonable assumption for winds near their launch 
point from a rotating accretion disk (Section 1). However, 
beyond the launch radius, if there are only radial forces, 
the conservation of angular momentum will lead to a ~l/r 
decline in the outflow transverse speeds. Magnetic fields 
threa ded through t he disk and applying a torque to the wind 
(e.g.. lEverettll2005h could counteract this decline. However, 
overall, we expect the actual transverse speeds to be less 
than or equal to the local Keplerian speed and, therefore, the 
distances derived from the Keplerian assumption are upper 
limits. (See below for more discussion.) These distances are 
shown as the left ordinate in Fig. 1151 This figure highlights 
the importance of obtaining measurements of BAL variabil- 
ity over shorter time-scales for constraining the location of 
the outflow. For example, variability times less than ~0.2 
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yr (knife-edge) or 0.6 yr (crossing disk) indicate distances 
smaller than the nominal radius of the C IV broad emission 
line region (BELR). 

The shortest time-scales over which we detect vari- 
ability is 0.02-0.03 yr (8 and 10 days) in 1246-0542 and 
0842+3431, with AA of 0.05-0.06. The transverse velocity 
of the outflow in these cases is 17 000 to 84 000 km s _1 
for 1246-0542 and 18 000 to 71 000 km s" 1 for 0842+3431, 
with the lower and upper limits given by the knife-edge and 
crossing disk models, respectively. These transverse veloci- 
ties are surprisingly high, from ~0.7 to 5 times the radial 
flow velocity, leading to derived distances of only ~0.001 to 
0.02 pc for both quasars, using the average SMBH mass for 
the quasar sample. For comparison, the estimated radius of 
the UV continuu m source is 0.004 pc (see above). Similarly, 
iMoe et al l i|2009l ) estimate a crossing speed of 18 000 km 
s _1 and a distance from the central SMBH <~0.1 pc, based 
on the variability in a single quasar. 

These results provide much smaller distance constraints 
than the changing ionization scenario discussed above (Sec- 
tion [S23J- There is indirect evidence for the crossing cloud 
scenario applying to 1246—0542 because its variations at 
— 17 200 km s _1 occurred in what appears to be an optically 
thick trough (Figs [2] and [3J). This interval lies in the bottom 
of the C IV BAL trough, at the same velocity as the bottom 
of the Si IV BAL trough. In Paper 2, we described how the 
Si IV optical depth should be ~3.4 times less than C IV, given 
the lower abunda nce of Si compared to C when assuming 
solar abundances (Hamann et al.lll997l ; iHamann fc Ferlandl 
ll999l ; lAsplund et all 2009). Furthermore, if the ionization is 
at least as high as that needed for a maximum C iv/C ratio, 
then the optical depth in Siiv should be >8 times smaller 
than C iv l|Hamann et aljfeoill ) . The apparent optical depth 
of Si IV in this interval is ~0.5, indicating that the optical 
depth of C IV is at least ~1.7— 4. This is greater than the ap- 
parent Civ optical depth of ~1.2, which suggests that the 
C iv line is saturated and would not likely be affected by 
modest changes in the ionization of the gas. 

It is interesting to note that the crossing speeds im- 
plied by the shortest variability times are similar to or sev- 
eral times larger than the observed radial outflow speeds. 
In models of thes e flows driven by radiation pressure 
jMurrav et al.lll995h . crossing speeds near or above the ra- 
dial speed will occur only if the flow we measure is near its 
launch radius, e.g., where the gas is still roughly in corota- 
tion with t he disk and the radial acc eleration is not yet com- 
plete (also. [Murray fc ChiangHl998T l. In models with signif- 
icant magnetocentrifugal forces, there can be large compo- 
nents of both azimuthal and vertical velocity (perpendicula r 
to the disk plane) far from the launch radius 
However, in luminous quasars with high accretion rates (i.e., 
high L/Ledd) , radiative (and thus radia l) forces are expected 
to dominate |Everettl2007i ; |Proga|[2003l 'l. Therefore, the most 
likely explanation for the large crossing speeds in our study 
is that the measured flows are, indeed, at very small dis- 
tances near their launch radius. 

One final note on the distances that we calculate for 
BAL outflows is that previous studies have foun d evidence 
for ab sorbers located beyond the BEL region. iTurnshekl 
(1988) describe cases where a BAL overlaps a BEL in the 
spectra, and the absorption goes deeper than can be ex- 
plained by the BAL absorbing just the continuum under- 



neath the BEL. This indicates that, at least for these par- 
ticular cases, the BAL region occults the BEL region. In our 
data, the Nv BAL in 0842+3431 might be deep enough to 
absorb both the continuum and some of the red side of the 
broad Lya emission line. This would place the BAL gas far- 
ther out than the BEL region, at radii r > 0.15 pc, and ap- 
parently contradict our estimates above for r ~ 0.001 — 0.02 
pc based on the variability and crossing speeds. However, 
these assertions need further study because i) very few 
quasars show evidence for the BAL gas occulting the BELs, 
and ii) in the particular case of 0842+3431, the extension of 
the red wing in the Lya emission profile is highly uncertain 
and the shape and depth of the observed N v trough might 
be significantly contaminated by absorption in the Lya for- 
est. 
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